//
//  InitialValue.cpp
//  NSESolver
//
//  Created by Patrick Zulian on 3/29/11.
//  Copyright 2011 ICS - USI - Lugano. All rights reserved.
//


#include "InitialValue.h"
#include "Types.h"
#include "CFDProblem.h"

namespace nsesolver {
    
    void InitToOnes(Matrixd &m)
    {
        m.allSet(1);
    }
    
    void InitToZeros(Matrixd &m)
    {
        m.allSet(0);
    }
    
    
    InitialValue::InitialValue(InitM u0, InitM v0)
    : u0(u0), v0(v0)
    { }

    
    InitialValue::~InitialValue()
    {
        //TODO
    }
    
    
    
    void InitialValue::setInitialValues(CFDProblem &problem)
    {
        Matrixd &u(problem.U());
        Matrixd &v(problem.V());
        
        u0(u);
        v0(v);
    }
}
